nd scattering lengths from a quark-model based NN interaction 
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^ ' Abstract 

We calculate the doublet and quartet neutron-deuteron scattering lengths 
using a nonlocal nucleon-nucleon interaction fully derived from quark-quark 
interactions. We use as input the NN 1 Sq and 3 S±- 3 Di partial waves. Our 
fN| \ result for the quartet scattering length agrees well with the experimental value 

while the result for the doublet scattering length does not. However, if we 
take the result for the doublet scattering length together with the one for the 
triton binding energy they agree well with the so-called Phillips line. 
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In a previous paper [1] the properties of the bound state of three nucleons, the triton, 
were studied using a nonlocal nucleon-nucleon (NN) interaction obtained from the chiral 
quark model [2]. In this Brief Report we would like to complete that study by calculating 
within the same framework the neutron-deuteron doublet and quartet scattering lengths. 

The nonlocal nucleon-nucleon interaction, obtained from the chiral quark model, is based 
on a Lippmann-Schwinger formulation of the resonating group method in momentum space. 
A detailed description of the method and a summary of the parameters of the model is given 
in Ref. [1]. In this reference it was shown that this NN interaction describes very well the 
properties of the deuteron as well as the 1 S , and 3 S 1 - 3 D 1 phase shifts. Besides, it was found 
that the nonlocal NN interaction obtained from the chiral quark model predicts a triton 
binding energy of 7.72 MeV, comparable to the predictions by conventional meson-exchange 
models. Thus, it is important to see if that result extends also to the neutron-deuteron 
scattering lengths. 

We will write down the integral equations that determine the nd scattering lengths using 
the partial-wave basis states 

I PiQi, Pi) =\ PiQi, ZiSijiUXiJi), (1) 

where if <7j and r, stand for the spin and isospin of particle i then £j, Sj, ji, ii, \, and Jj 
are the orbital angular momentum, spin, total angular momentum, and isospin of the pair 
jk while Aj is the orbital angular momentum between particle i and the pair jk and Jj is 
the result of coupling A, and <7j. The conserved quantum numbers are J, the total angular 
momentum and I, the total isospin. 

The Faddeev equations that determine the nd scattering lengths in the special case when 
one considers only S wave configurations were written down in Ref. [3]. The generalization 
of these equations to the case when one includes arbitrary orbital angular momenta are 

(PiQhPi I T/ 1 | O ) = 25 Stl 5 jil 5 uo Sx t o\s{q i )GQ 1 {E;p i q i )(l) ei (Pi) 

Qi 

roo rl 3^? 

+ S / 1j d ^ J dcos9 tWi'diU (Pi,Pi, E - 4]^) 

e 'iPi 

x G (E; I'j'l,)!);]''; ,('!;. qj, cos#) (p^-; Pj \ Tf 1 \ <f> ), (2) 

where t^i'.s^iXPiiP'i'i e ) are the nucleon-nucleon t-matrices, M is the mass of the nucleon, 
and E = —B d where B d is the binding energy of the deuteron. is the deuteron wave 

function with orbital angular momentum £i, and 

G (E; Pt q t ) = E _ p 2 /M ^ 3q?/4M + te , (3) 
Pi = ^/q] + qf/^ + q l q j cos9, (4) 
Pj = \/<li + q] / 4 + q>q jCOsO, (5) 
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Pi = {V'iSijiiiXiJi}. 



(6) 



D i ^.J I (q i ,qj,cos9) are the angular momentum-spin-isospin recoupling coefficients which are 
defined by Eqs. (21)-(26) of Ref. [4]. 

After obtaining the solution of Eqs. (2) the neutron-deuteron scattering lengths are 
calculated as 

2Mvr 



aj = — =— 2^ 2^ ^i^i^-oOAj-o / % dQi 



Pj Pi 



x 0,.(ft)Z^\ (0, ft, 0)(ft/2, ft; Pl | Tp | O ). (7) 

Since in Ref. [1] the triton binding energy was calculated using as input the NN 1 S and 
3 S\- 3 Di partial waves we will use the same prescription in our calculation of the neutron- 
deuteron scattering lengths. This implies a five-channel Faddeev calculation for the doublet 
scattering length ai/ 2 and a seven-channel calculation for the quartet scattering length a 3 / 2 . 
We give these channels in Table I. 

Our method [4] to solve the Faddeev equations consists in transforming them from being 
integral equations in two continuous variables into integral equations in just one continuous 
variable. This is achieved by expanding the two-body t— matrices in terms of Legendre 
polynomials as 

U(pi,tie) = T, P n(^rr(e)PM), (8) 

nr 

where P n and P r are Legendre polynomials, 

Pi + o 



and pi and p\ are the initial and final relative momenta of the pair jk while b is a scale 
parameter on which the results do not depend. We found that using 6 = 3 fm _1 leads to 
very stable results while for the expansion (8) we found convergence with twelve Legendre 
polynomials, i.e., < n < 11. 

We have calculated also the triton binding energy to check that our calculation reproduces 
the value obtained in Ref. [1]. We give in Table II our results for the triton binding energy 
and the two neutron-deuteron scattering lengths as well as the corresponding experimental 
results [5]. As pointed out in Ref. [1], the result B = 7.72 MeV for the triton binding energy 
predicted by the chiral quark model is comparable to the values obtained by conventional 
meson-exchange models like Nijmegen or Bonn [6] since the theoretical value differs by less 
than 1 MeV from the experimental result. The situation in the case of the doublet scattering 
length ai/2 appears to be somewhat worse since the theoretical value 1.13 fm is almost a 
factor of two larger than the experimental result. However, as we will see next, that is not 
the case. As it is well-known [7], the results obtained from a given theoretical model for the 
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triton binding energy B and the doublet scattering length au2 are strongly correlated. The 
results obtained from different models follow what is known as the Phillips line which is a 
straight line relating B versus a\/2- For example, in Ref. [8] a five-channel calculation similar 
to ours was performed using three local potentials [9-11]; therefore, we made a minimum- 
square fit of their results for B and a±/ 2 to obtain the Phillips line ai/ 2 = 6.352—0.6775 which 
for B = 7.72 MeV gives ai/2 = 1.13 fm in agreement with our result. The most complete 
calculations of B and a±/2 have been performed in Ref. [12] for a variety of modern nucleon- 
nucleon force models [13-15] and where they have included the higher angular momentum 
two-body channels as well as three-body forces [16-19]. Using the values of B and ai/ 2 for the 
48 different models considered in Ref. [12] we obtained the Phillips line 01/2 = 7.028 — 0.7565 
which for B = 7.72 MeV will give a doublet scattering length ai/2 = 1.19 fm which is also 
quite close to the result of our calculation. Our result for the quartet scattering length 
a 3/2 = 6.40 fm agrees well with experiment and with Refs. [8,12]. 

The reason why our calculation agrees well with experiment in the case of the J = 3/2 
channel but not in the case of the J = 1/2 channel is that the former is determined by 
the pure S'-wave configuration £ = A = while the latter has important contributions from 
higher angular momentum two-body channels and from three-body forces, both of which 
are lacking into our model. This is very similar to the situation encountered in effective 
field theory where the J = 3/2 channel is very well explained by S'-wave models without 
three-body forces [20,21] while for the J = 1/2 channel the theory can produce sensible 
results only when a three-body force is included [22,23]. 

In summary, we conclude that the nonlocal NN interaction obtained from the chiral 
quark model gives results for the neutron-deuteron scattering lengths which are comparable 
to those obtained from conventional meson-exchange models. 

This work has been partially funded by Ministerio de Education y Ciencia under Contract 
No. FPA2007-65748 and by COFAA-IPN (Mexico). 
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TABLES 



TABLE I. Three-body channels that contribute to a given NNN state with total isospin / and 
total angular momentum J. 
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TABLE II. Triton binding energy B (in MeV) and scattering lengths and a 3 / 2 (in fm) as 
compared with the corresponding experimental values. 



Quantity 


Theory 


Experiment 


B 


7.72 


8.48 


ai/2 


1.13 


0.65±0.04 


a 3/2 


6.40 


6.35±0.02 
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